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Comment on: Thermostatistics of Overdamped 
Motion of Interacting Particles 

In a recent Letter [l[ Andrade et al. argued that over- 
damped particles in contact with a reservoir at T = 
"obey Tsallis statistics with entropic index v — 2". To 
justify this, Andrade et al. performed a simulation of 
classical particles in two dimensions, interacting through 
a pair potential V(r) = gG(xi,X2), where G(xi,X2) = 
qKo(\x.i — x 2 |/A), and K is a modified Bessel function of 
order zero. This potential models a vortex-vortex inter- 
action in a type II superconductor, where q is the vortex 
strength and A is the effective London penetration length. 
The particles were confined in a Id parabolic trap with 
a potential W(x) = ^x 2 , while a periodic boundary con- 
dition was used in the y direction (periodicity L y ). The 
simulation was performed using overdamped molecular 
dynamics (MD) with the thermostat set at T = 0. The 
resulting stationary density profile was then fitted to a 
distribution that maximizes the Tsallis entropy of index 
v = 2. 

The idea that non-extensive Tsallis statistics is rele- 
vant for particles at T = is very surprising and, if 
correct, would certainly justify publication in Physical 
Review Letters. However, prior to discarding the stan- 
dard statistical mechanics one should see what it has to 
say on this matter. 

Within traditional thermodynamics, if a classical sys- 
tem is placed in contact with a temperature reservoir at 
T = it will loose all its kinetic energy and collapse to 
the ground state — the minimum of the potential en- 
ergy. In the ground state the net force on each particle 
vanishes. We will now show that this is precisely what 
happens for the system studied by Andrade et al. 

From now on we will measure all lengths in units of A. 
We will work in thermodynamic limit N — > oo, at fixed 
total vertex strength, q 2 N — 1. We note first that the 
potential produced by a particle located at xi satisfies 



V 2 G(x, xi) - G(x, Xl ) = -27rg(5(x - Xl ; 



(1) 



Taking into account the periodicity in the y-direction, 
this equation can be solved exactly to yield 



G(x;x 1 ) = 7 £ 



,-~/m\X-Xl\ 



In 



(2) 



where 7 m = y 1 + Air 2 m 2 / L 2 . Furthermore, symmetry 
requires that at equilibrium the total potential inside the 
system is a function of ^-coordinate only, tp(x). The force 
balance on each particle then reduces to q<p'(x) — —ax. 
In view of eq. ([1]), f{x) must also satisfy the inhomo- 
geneous Helmholtz equation, <p"(x) — ip(x) = —2irqp(x), 
from which we conclude that the particle density p(x) is 
a restricted parabolic function, 



p{x) 



Airq' 



Q(x 2 



(3) 



where O is the Heaviside step function. Note that the 
density does not go to zero smoothly; instead it extends 
up to ±xo, after which it drops discontinuously to zero. 
The value of x m is determined by the density normal- 
ization, N = J p(x)alxdy, x m = \J a "^ Xg ~ and the 
extent of the density distribution by the force balance, 
qip'(xo) = —axo. Using the Green function ([2]), the force 
balance on a particle at x — xq requires 



1 



(xl l -x 2 )e- ( - x °-^dx=x , 



(4) 



which reduces to, 



-2.r 



(2 + 2x Q + x 2 - x 2 m ) =4:ro. 
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FIG. 1: The vertex strength density distribution p q (x) 
q 2 p(x) for a 



10 and L v 



20. Circles are the result 



of our MD simulation with iV = 1000 particles, while the 
solid curve is the prediction of the present theory. There are 
no adjustable parameters. 

This constitutes the exact solution for the particle dis- 
tribution at T = 0. In Fig. 1 we compare it with the 
MD simulation. As expected, an excellent agreement is 
found between the theory and the simulation, without 
any fitting parameters. We conclude, therefore, that the 
density distribution of particles in contact with a reser- 
voir at T — has nothing to do with the Tsallis statis- 
tics, and everything to do with Newton's Second Law. Of 
course for finite temperature, the density distribution will 
be described by the usual Maxwell-Boltzmann statistical 
mechanics. 
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